library(openxlsx)
library(stats)

dat_pca<-read.xlsx("./data.xlsx")
pcaresults<-prcomp(dat_pca,center = T,scale. = T)

fviz_pca_var(pcaresults,col.var="black",col.circle = "black",repel = T,axes = c(1,2),geom.var = c("arrow","text"))+theme_bw()+theme(text=element_text(family="Times New Roman",face="bold",size=18))+xlab("PC 1 (34.1%)")+ylab("PC 2 (25.1%)")
fviz_pca_var(pcaresults,col.var="black",col.circle = "black",repel = T,axes = c(1,3),geom.var = c("arrow","text"))+theme_bw()+theme(text=element_text(family="Times New Roman",face="bold",size=18))+xlab("PC 1 (34.1%)")+ylab("PC 3 (20.6%)")